#include "Interpolation.h"

int main()
{
    vector<double> D;
    vector<double> Sp1;
    vector<double> Sp2;
    D = {0, 6, 10, 13, 17, 20, 28};
    Sp1 = {6.67, 17.3, 42.7, 37.3, 30.1, 29.3, 28.7};
    Sp2 = {6.67, 16.1, 18.9, 15.0, 10.6, 9.44, 8.89};
    Newton N1(D, Sp1);
    Newton N2(D, Sp2);

    cout << "D = " << "";
    for (int i = 0; i < D.size(); i++)
    {
	cout << D[i] << " ";
    }
    cout << endl;
    cout << endl;
    cout << "Sp1 = " << "";
    for (int i = 0; i < Sp1.size(); i++)
    {
	cout << Sp1[i] << " ";
    }
    cout << endl;
    N1.printF();

    cout << "Sp2 = " << "";
    for (int i = 0; i < Sp2.size(); i++)
    {
	cout << Sp2[i] << " ";
    }
    cout << endl;
    N2.printF();

    for (int i = 0; i <= 50; i++)
    {
	cout << "-" << "";
    }
    cout << endl;
    return 0;
}
